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Abstract 

The Crigler-Najjar Syndrome Type I (CNSI) is a rare genetic disorder caused by mutations in the Ugtlal gene. It is 
characterized by unconjugated hyperbilirubinemia that may result in severe neurologic damage and death if untreated. To 
date, liver transplantation is the only curative treatment. With the aim of generating mutant cell lines of the Ugtl gene, we 
utilized the TALEN technology to introduce site-specific mutations in Ugtl exon 4. We report a fast and efficient method to 
perform gene knockout in tissue culture cells, based on the use of TALEN pairs targeting restriction enzyme (RE) sites in the 
region of interest. This strategy overcame the presence of allele-specific single nucleotide polymorphisms (SNPs) and 
pseudogenes, conditions that limit INDELs' detection by Surveyor. We obtained liver-derived murine N-Muli cell clones 
having INDELs with efficiency close to 40%, depending on the TALEN pair and RE target site. Sequencing of the target locus 
and WB analysis of the isolated cell clones showed a high proportion of biallelic mutations in cells treated with the most 
efficient TALEN pair. Ugt glucuronidation activity was reduced basal levels in the biallelic mutant clones. These mutant liver- 
derived cell lines could be a very useful tool to study biochemical aspects of Ugtl enzyme activity in a more natural context, 
such as substrate specificity, requirement of specific co-factors, the study of inhibitors and other pharmacological aspects, 
and to correlate enzyme activity to the presence of specific mutations in the gene, by adding back to the mutant cell clones 
specific variants of the Ugtl gene. In addition, since genome editing has recently emerged as a potential therapeutic 
approach to cure genetic diseases, the definition of the most efficient TALEN pair could be an important step towards 
setting up a platform to perform genome editing in CNSI. 
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Introduction 

UDP-glucuronosyl transferases (UGTs) catalyze glucuronida- 
tion of a great variety of compounds. They are classified into two 
subfamilies, UGT1 and UGT2 [1]. Each subfamily is composed of 
several different isoforms, expressed in a tissue-specific manner 
and having substrate specificity [2] . The most important member 
of the UGT 1 subfamily is UGT 1 a 1 , which is mainly expressed in 
the fiver, and is the only enzyme able to conjugate bilirubin. 
Deficiency of this isoform results in Crigler-Najjar Syndrome Type 
I (CNSI), a rare genetic disorder characterized by unconjugated 
hyperbilirubinemia [3,4]. Persistent unconjugated hyperbilirubin- 
emia may cause severe neurologic damage and result in death by 
kernicterus if untreated. To date, liver transplantation is the only 
curative treatment. 

The study of the mechanisms of disease and the effects of 
pharmacological therapies rely in the use of specific cellular and 
animal models containing mutant versions or variants of the 
causative gene. The absence of rapid and efficient tools to modify 
the genome limited the availability of cellular models. However, in 
the last few years the development of a new generation of 
nucleases having a very high sequence-specificity allowed the rapid 
and efficient modification of the loci of interest [5-8]. These 



nucleases are known as Zinc Finger nucleases (ZFN), Transcrip- 
tion Activator Like Effector Nucleases (TALEN) and Clustered 
Regularly Interspaced Short Palindromic Repeats (CRISPR) 
together with the CRISPR associated nuclease 9 (Cas9), also 
denominated RNA Guided Endonucleases (RGEN). They gener- 
ate a precise single- or double- strand break at a specific locus 
[6,9]. Among the mentioned nucleases, TALENs integrate two 
important features: the ease of construction that ZFN lack, and 
higher specificity than RGENs have produced so far. They consist 
of two domains: a customizable DNA binding domain and the 
catalytic domain of Fokl endonuclease [7,9]. TALENs work as 
heterodimers in which each monomer binds 15-20 bp of DNA 
that flanks a 15-24 bp spacer region. Once both monomers bind 
to the target sequence, the Fokl endonuclease domains introduce a 
double strand break (DSB). When homologous donor DNA is 
present, DSBs are repaired by Homology-directed repair (HDR), a 
pathway that corrects the damage and is used for "gene editing". 
Alternatively, when no donor DNA is present, the damage is 
repaired by non-homologous end-joining (NHEJ), which results in 
insertions or deletions (INDELs) in the target site [6,7]. 
Consequendy, when the target site is in the protein coding 
sequence, two-thirds of the repairs will cause a frameshift, which 
may lead to translational termination and loss of function [6,7]. 
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In the present work, we used the TALEN technology to 
introduce mutations in the Ugtl gene of mouse liver-derived N- 
Muli cells. Our strategy was based in the inactivation of the gene 
by targeting restriction enzyme sites present in the exon 4 of the 
Ugtl. We successfully generated mouse liver cells containing 
mono- and bi-allelic mutations in the exon 4 of the Ugtl gene. 

Materials and Methods 

Generation of TALENs pairs 

To design the different TALEN pairs we used the online TAL 
Effector Nucleotide Targeter (TALENT) software [10]. To 
generate the TALEN plasmids we used the "Golden Gate 
TALEN and TAL effector Kit 2.0 (Addgene, cat. 
#1000000024), as previously described [11]. The generated 
TALENs are based on the Goldy TALEN scaffold, that contains 
a deletion of 152 amino acids (aa) at the N-terminal domain and a 
63 aa C-terminal residue [12]. This scaffold was previously shown 
to be the most efficient one [9]. The target sequences of the 
TALENs are indicated in Figure S3. 

Surveyor assay 

To set up Surveyor Assay genomic DNA (gDNA) was extracted 
from wild type (WT) and UGT1A1 KO animals by overnight 
proteinase K digestion of tail biopsies. Plasmids containing WT 
and mutant (single point mutation) sequence were used as control. 
Several primer pairs were tested for both gDNA and plasmid 
DNA, with the highest yield and a single PCR product amplified 
with the LSur 1 F and LSur 1 R primer pair. For the Surveyor Assay 
the DNA was then PCR-amplified with LSurlF and LSurlR 
primers pair (see Table SI) under the following conditions: 95°C 
5 min, then 94°C 30 s, 70°C 30 s, 72°C 30 s for 35 cycles, with 
3 min of final extension step at 72°C. After PCR the amplicons 
were treated with Surveyor (Surveyor mutation detection kit, 
Transgenomic) as per manufacturer's instructions. Products were 
resolved by agarose gel electrophoresis (2% agarose gel). 

Generation of the luciferase reporter vector 

The pGL3-Linker vector was obtained from the pGL3-control 
vector by inserting the 582 bp long Ncol-Bcll fragment (purified 
from the pPGL-3 control vector) and a 3 1 bp linker (containing the 
BclI-EcoRI-EcoRV-BamHI-Pstl-Narl restriction enzyme sites) 
into a pGL-3 vector (having the BamHI site inactivated) previously 
cut with Ncol and Narl restriction enzymes. The final construct 
was verified by restriction enzyme digestion and sequencing. 

To generate the pGL3-mUGTlal WT vector, the 706 bp long 
BamHI-BamHI fragment of the mUgtlal-Ex4 WT, (previously 
cloned and sequenced in our lab), was cloned into the BamHI site 
of the pGL3-linker plasmid. 

Transient cell transfections and dual luciferase assay 

24 h before transfection HEK293 cells were plated in 24 wells/ 
plate in complete DMEM (10%FBS and antibiotic). At the day of 
transfection the cells were at 50-70% of confluence. The cells were 
co-transfected with 40 ng of Renilla (phRG-TK, Clontech), 1.6 |xg 
of pGL3-Linker (as negative control) or pGL3-mUGTlal-ex4 
plasmids and increasing amounts of the TALEN plasmids pair 
(400 ng, 800 ng and 1.6 u.g) using Lipofectamine 2000 as 
described by the manufacturer (InVitrogen). The day after, 
500 ul of complete DMEM (containing 20% FBS) were added. 
72 h after transfection cells were harvested and TALEN cutting 
efficiency was indirectly determined by quantifying luciferase 
activity with the Dual Luciferase assay (Promega). Briefly, cells 
were washed and lysed with 100 u.1 of PLB 1 x (Passive Lysis Buffer 



-PROMEGA), incubated 15 min at RT and centrifuged at top 
speed for 1 min. Firefly Luciferase activity was first measured by 
adding 30 ju.1 of LAR II (Luciferase Assay Reagent II) to 30 ul cell 
lysate; then 30 ul of Renilla luciferase substrate (Stop & GLO) 
were added and renilla luciferase activity was determined with a 
single-sample 20/20" luminometer from Turner Biosystem. 

Generation of mutant stable cell clones by TALEN 
transfection 

N-Muli cells (ATCC CRL-1638) were cultured in D-MEM with 
10% FCS. Due to the low transfection efficiency of N-Muli cells 
(15-20%), we performed two successive rounds of transfection: a) 
2xl0 6 N-MuLi cells were transfected with 2 ng of p-EGFP-C2- 
plasmid (for G418 selection, Gentech) and 3 ng of both TALEN 
right and left arms cloned in the pC-Goldy plasmid (Mnl, Nco, Nla 
or BstXI), by using the Effectene reagent (Qiagen) and following 
the manufacturer's instructions. Cells were cultured at 37°C for 
24 h (with 10% FBS) and 72 h later G418 (250 u.g/u.1) was added; 
b) When cells reached 70-80% of confluence they were transfected 
as indicated above, but with 2 ug of pMCS-AAT-PB:PGKpuro- 
deltaK (Welcome Trust Sanger Institute) for puromycin selection, 
instead of the p-GFP-plasmid. 72 h after transfection Puromycin 
at 2 ng/ul was added and single resistant clones were selected, 
picked and individually cultured. 

Genomic DNA preparation, PCR of genomic DNA, 
digestion with restriction enzymes to analyze INDELs, 
and DNA sequencing. 

Genomic DNA preparation was performed as previously 
described [13]. Briefly, cells were lysed in 50 mM TRIS pH 7.5, 
100 mM EDTA, 0.5%SDS and 200 ug/ml of protease K for 1 
hour at 37°, centrifuged for 10 min at maximum speed, genomic 
DNA precipitated with 2-propanol and resuspended in TE after 
centrifugation and washing. 

Genomic PCR was performed from 10-30 ng of DNA, by using 
the primers indicated in Table SI. The PCR product was then 
digested with the restriction enzymes specific for each TALEN 
pair (Nlal, Nco, BstXI and Mnll), and the RE-resistant fraction 
was purified, cloned into pUC19 and sequenced. 

Western blot analysis of positives clones 

Protein cell extracts were prepared as already described [13]. 
Briefly, cells were lysed in 50 mM Tris pH7.5, 150 mM NaCl, 1 % 
NP40 (containing protease inhibitors, Complete-Roche), centri- 
fuged to clarify debris, and protein concentration was quantified 
by Bradford (Biorad). The same amounts of protein extracts 
(100 ug for cell clones, and 30 for liver) were loaded on 10% SDS- 
PAGE, blotted onto nitrocellulose, incubated with the primary 
antibody anti-mouse UGTlal (Santa Cruz, cat# sc-25847, 
1:400), incubated with a secondary antibody (anti-rabbit HRP- 
conjugated - DAKO) and developed by ECL (Pierce) as described 
[13]. 

Ugt enzyme activity assay 

The UGTlal activity was determined by UGT-GloTM assay 
kit in the microsomal samples according to manufacturer's 
instructions (Promega, Madison, WI). The rationale of the assay 
resides in the use of a UGT Multienzyme Substrate, which reacts 
with the luciferin detection reagent to give light that can be 
quantitated by a luminometer. When the UGT substrate is 
glucuronidated it will no longer react with the luciferin detection 
reagent to give light. To determine enzyme activity, two reactions 
are set up per sample, one with the co-substrate UDPGA and the 
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other without. Only the reaction with the UDPGA will produce 
glucuronidated substrate if UGTs are present and active. Total 
UGT activity is then quantified by the difference in emitted light 
between the two reactions. To determine UGT1 activity, we 
prepared microsomes from wt N-Muli, Mnl20, Nco20 and Nco25 
cells. We used liver microsomes prepared from Ugtl KO mice and 
WT mice as controls [13], in addition to the microsomal 
preparations (positive control) of the UGT-GloTM assay kit 
(Promega Madison, WI). Microsomes (10 ug from cell clones, 1 u,g 
from mouse liver and 4 ug from positive control) were incubated 
with 16 mM UDPGA, 20 U.M multi enzyme substrate and 
0.35 mg/ml digitonin (Sigma) at 37°C for 150 min. After which 
40 m.1 of reconstituted luciferin detection reagent containing D- 
Cysteine was added and the luminescent signal was allowed to 
stabilize for 20 min at room temperature. Luminescence was read 
using a multiplate reader (Perkin Elmer Envision Plate Reader, 
Walthan, MA, USA). Data was expressed as percentage of the 
luminescent substrate consumed. 

Results 

Targeting restriction enzyme sites with TALENs to avoid 
interference of strain-specific Ugtl SNPs and/or 
pseudogenes 

In order to determine TALEN efficacy in targeting the mUgtl 
gene, we initially performed the Surveyor assay using PCR 
amplification products obtained from plasmids containing the wt 
or mutated sequences of the Ugtl gene exon 4 and flanking 
introns. These DNA fragments, originating from 1 29Sv/J ES cells, 
had been previously cloned and mutated in our laboratory, and 
used to generate a mouse strain containing one base deletion in the 
exon 4 of the Ugtl gene [13]. To our expectations, the Surveyor 
assay resulted in two bands of the predicted size when wt and 
mutant PCR products were mixed in the presence of Surveyor 
(192 and 324 bp, Figure SI, lanes 2-4). Surprisingly, when we 
performed the Surveyor assay using the PCR products amplified 
from wt and mutant Ugtl mice genomic DNA, we obtained a 
laddering pattern of bands, regardless of the source of DNA used 
for the reaction, indicating the presence of heteroduplexes (Figure 
SI, lanes 5-7). A similar result was obtained when we tested 
genomic DNA from N-Muli hepatoma cells (data not shown), 
which originates from the C57B1/6 mouse strain. 

Since mutant mice had been backcrossed with C57B1/6 mice 
for 9 generations to obtain a homogeneous genetic background, 
but the mutation was generated in 129Sv/J ES cells, we suspected 
the presence of strain-specific variants and/ or pseudogenes that 
might produce mismatches and unexpected bands in the Surveyor 
reaction. To verify this possibility, we cloned the PCR product of 
wt 129Sv/J genomic DNA and sequenced nine (9) independent 
clones. As reported in Figure S2, we compared those clones with 
the available C57B1/6 GenBank reference sequence and we 
observed the presence of mismatches among some of the 129Sv/J 
clones with the reference C57B1/6 sequence and, unexpectedly, 
among the 129Sv/J-derived clones themselves. The sequences 
were categorized into two distinct groups according to the 
mismatches, as indicated in Figure S2. We speculate that the 
failure in the Surveyor assay could be due to the presence of strain- 
specific polymorphisms in the Ugtl exon 4, and/or to a Ugtlal 
pseudogene. The presence of a pseudogene of the Ugtl gene was 
previously suggested by data presented by Miles [14], Seldin 
[15,16], Lathman et al [17], MGI and NCBI (UDP glucuronosyl- 
transferase 1 family, polypeptide Al, related sequence 1, Ugtlal- 
rsl, Jackson Lab, http://www.informatics.jax.org/marker/ 
MGL98899; NCBI, http://www.ncbi.nlm.nih.gov/gene/ 



111282). However, no active link to the pseudogene sequence 
was yet available in the mentioned gene banks and we failed to 
identify the Ugtlal-rsl pseudogene in the currendy available gene 
banks (not shown). 

Therefore, to tackle the problems in setting up the Surveyor 
assay, we modified the strategy by targeting restriction enzyme 
(RE) sites present in the exon 4 of the mouse Ugtl gene. TALEN- 
induced double strand break will result in the disruption of the RE 
site by error-prone NHEJ repair. In this case, digestion of the PCR 
product with the RE will result in a fraction of the PCR product 
that will be RE-resistant. We designed and constructed four (4) 
TALEN pairs targeting restriction sites present in the exon 4 of the 
Ugtl gene (NlalV, Ncol, BstXI and Mnll, Figure 1). The selected 
restriction sites were located in the spacer region, between the two 
TALEN monomers. The sequences recognized by the different 
TALENs are indicated in Figure S3. To construct the TALEN 
plasmids we used the "Golden Gate TALEN and TAL effector Kit 
2.0 [11], based on the Goldy TALEN scaffold [12] 

Testing TALEN pairs activity with a luciferase reporter 
plasmid 

To determine the efficiency of the different pairs of TALENs, 
we constructed a reporter vector containing a duplicated region of 
the firefly luciferase cDNA, with the exon 4 of the Ugtl gene 
inserted in between the repeated luciferase regions (Figure 2A). 
Double strand breaks generated by TALENs in the exon 4 will 
favor the recombination events within the luciferase gene and, in 
turn, restore the correct primary sequence in the reporter vector. 
Efficiency of the different TALEN pairs will be then determined as 
luciferase activity. 

We first transfected N-Muli cells with the pGL3-(Linker) 
reporter vector without insert and verified the absence of luciferase 
activity (Figure 2B). Then, we co-transfected the reporter vector 
containing the Ugtl exon 4 [pGL3-UGTlal-(Ex4)-WT] together 
with the different TALEN pairs and determined the activity of the 
reporter plasmid. Transfection of all four TALEN pairs resulted in 
sustained luciferase activity, with the NcoI-TALEN being the most 
active one, followed by the Mnll-TALEN. The NlalV- and BstXI- 
TALEN pairs showed lower luciferase activity, that amounted to 
about half of the activity observed for the Mnll- and NcoI-TALEN 
pairs (Figures 2C and 2D). The renilla-luciferase vector was used 
to normalize for differences in transfection efficiency. 

Transient transfection of Ncol TALEN pair results in short 
deletions at the Ncol site 

We disrupted the Ugtl gene in N-Muli mouse hepatoma cells 
by transiently transfecting the most effective TALEN pair, the 
NcoI-TALEN. N-Muli cells were selected because Ugtl expression 
levels were high enough to be detected by Western blot (see 
below). Due to the low transfection efficiency of N-Muli cells 
(between 15 and 20%, data not shown), we used a vector 
expressing the GFP protein, together with the TALEN plasmids, 
to co-transfect the cells. GFP-positive cells were sorted by FACS, 
genomic DNA was prepared, the exon 4 region of the Ugt 1 gene 
was PCR amplified and the amplicon analyzed by RE digestion. 

We observed that the PCR product amplified from GFP- 
transfected N-Muli cells was completely digested by Ncol while, in 
the cells transfected with the Ncol TALENs, a fraction of the PCR 
product was resistant to Ncol digestion (estimated in ~8%, 
Figure 3 A, lane 5). To further confirm the inactivation of the Ncol 
site by TALENs, the Ncol-resistant fragment was purified from 
the agarose gel, PCR amplified with an internal pair of primers 
(nested PCR), and the PCR product was digested with Ncol. As a 
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Figure 1. Scheme of the WT exon 4 showing the position of the different restriction enzyme sites and corresponding TALEN pairs. 

The murine Ugtl exon 4 and flanking introns is shown in the figure. The relative position of the NlalV, Ncol, BstXI and Mnll restriction sites is 
indicated, as well as the left and right arm of each TALEN pair. Exonic and intronic sequences are indicated in yellow and white, respectively. 
doi:1 0.1 371 /journal.pone.01 0481 6.g001 



negative control, we purified the PCR product originating from 
cells transfected only with the GFP plasmid. When the PCR 
products were digested with Ncol, we observed complete digestion 
of the band derived from the GFP-transfected cells, while 80% of 
the nested-PCR band originating from the TALEN-transfected 
cells was resistant to Ncol digestion (Figure 3B), indicating that the 
band shown in Figure 3A, Lane 5 was resistant to Ncol digestion 
due to inactivation of the Ncol site. To characterize the INDELs 
generated in the N-Muli cells, we cloned and sequenced the PCR 
product of the Ncol resistant fraction (Figure 3A, Lane 5). We 
observed that the majority of the mutations were deletions of the 
target region ranging from 3 to 1 1 bases, with the exception of 
three clones containing single base mutations (Figure 3C). On the 
contrary, no insertions were observed. All deletions and mutations 
were in the spacer region and, as expected, in all cases inactivated 
the Ncol restriction enzyme site. 

Isolation and screening of Ugtl deficient clones 

To overcome the obstacles caused by the low transfection 
efficiency of N-Muli cells, we adopted a two-steps strategy based in 
the co-transfection of the TALEN plasmids with a second vector 
encoding for a selection marker (see Figure S4 for a flow diagram 
of the method). First, we co-transfected the cells with the pEGFP- 
C2 plasmid, which encodes for the Neomycin selection marker, 
together with the TALEN vectors. G418-resistant cells were 
selected for 12-13 days (when G418 cells reached about 70% 
confluence). Second, G4 1 8-resistant cells were then transfected 
with the same TALEN plasmids and a second plasmid encoding 
for another selection marker (Puromycin gene). Transfected cells 
were grown in the presence of both antibiotics (G418 and 
puromycin) and individual clones were isolated and amplified. 

We tested the four TALEN pairs mentioned above (Ncol, Mnll, 
BstXI and NlalV) in N-Muli cells and analyzed 30, 30, 12 and 10 
G418/puromycin resistant clones of each transfection, respectively 



(Table 1). As an example, the analysis of two Ncol and two Mnll 
clones are shown in Figure 4. The region of interest was a 680 bp- 
long fragment and contained one Mnll site located in exon 4 of the 
Ugtl gene, which was targeted by the Mnll-TALEN, and two 
other Mnll sites located in the 3' flanking intron (Figure 4A). 
Complete digestion was observed when the amplicon obtained 
from untransfected N-Muli cells was treated with Mnll, while a 
fraction of the PCR product from Clone Mnl4 was resistant to 
Mnll digestion (Figure 4C) suggesting the presence of a mono- 
allelic mutation. In the case of Clone Mnl20, we obtained a shorter 
PCR product in addition to the expected PCR product, suggesting 
the presence of a larger deletion in one of the alleles. Both bands 
were digested by Mnll at the site present in the intron at the 3' end 
of the amplicon, but were not cut at the Mnll target site located at 
the 5' splicing site of Exon 4. These results suggest the presence of 
a large deletion in one of the alleles, and the inactivation of the 
Mnll site in the other allele. 

In the case of the Ncol RE approach, only one Ncol site was 
present in the 563 bp-long PCR fragment containing Ugtl exon 4 
(Figure 4B). Digestion of the wt fragment resulted in two bands of 
295 bp and 268 bp (Figure 4D). Digestion of the PCR product of 
Clone Nco25A showed the presence of a Ncol-resistant fraction of 
the amplicon, suggesting the presence of a mono-allelic mutation. 
Similarly to that observed for the Mnl20 clone, we observed two 
PCR products for the Nco20 clone. Both bands were resistant to 
Ncol digestion, suggesting the presence of a large deletion in one 
of the alleles, and the inactivation of the Ncol site in the other 
allele. As previously observed in the Mnl20 clone, we observed 
faint bands corresponding to the digestion products of the wt 
amplicon. 

Molecular analysis of Ugtl -deficient clones 

To characterize the genetic defect generated by TALENs we 
amplified, cloned and sequenced the Ugtl exon 4 and flanking 
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Figure 2. Evaluation of TALEN activity using a Luciferase-reporter plasmid. A) Scheme of the TALEN reporter plasmids and strategy. The 
Ugt1 exon 4 and flanking introns were introduced into the polylinker of the pGL3 (Linker) plasmid. This plasmid contains a 548 bp-repeated region of 
the luciferase cDNA. DSB induced by TALENs stimulates homologous recombination between the repeated luciferase regions and results in the 
recovery of luciferase enzyme activity. B) The reporter plasmid without TALENs showed no luciferase activity. 0.4 jig, 0.8 |ig and 1.6 |ig of 
pGL3-control and pGL3 (Linker) plasmids (lanes 1 -3 and 4-6, respectively) were co-transfected with the 40 ng of renilla-luc plasmid into HEK293 cells. 
The pGL3 (Linker) plasmids showed no luciferase activity (lanes 4-6). C) Scheme of the reporter plasmid used in the TALEN assays. The figure 
shows the schemes of the pGL3-control, pGL3 (Linker) and pGL3-Ugt1(Ex4)WT plasmids. D) The Ncol TALEN pair shows the highest activity. 
The Mnll, Ncol, NlalV and BstXI TALEN pairs were co-transfected into HEK293 cells together with the reporter pGL3-Ugt1 (Ex4)WT and renilla luc 
plasmids. Increasing amounts of the TALEN plasmids were co-transfected (200 ng for lanes 2, 6, 10 and 14; 400 ng for lanes 3, 7, 1 1, and 15; 800 ng 
for lanes 4, 8, 12, and 16). Activity was normalized to lanes 1, 5, 9, and 13, which contain pGL3 (linker) co-transfected with the higher amount of the 
respective TALEN pair. 
doi:1 0.1 371 /journal.pone.01 0481 6.g002 



introns from the isolated cell clones. Both in the Mnll and Ncol 
clones we observed the presence of deletions, most of them ranging 
from 2 to 1 1 bases (Figure 5A). However, we found larger 
deletions of 353 bp, 64 bp, 42 bp, 147 bp and 107 bp in the 
Mnl20, Mnl7, Ncol 5 (both alleles) and Nco20 clones, respectively, 
confirming the results obtained by PCR of genomic DNA. We 
found bi-allelic mutations in the Nco8, Ncol 5, Ncol 9, Nco20 and 
Nco27 clones (Figure 5A). 

Next, we determined Ugtl protein levels in the selected cell 
clones. We prepared protein extracts from the isolated clones and 
performed a WB analysis to determine the presence of the Ugt 1 
protein. We observed a significant reduction in protein levels in 3 
Mnll clones (Figure 5B, Mnl4, Mnll 5 and Mnl29, while no Ugtl 
band was detected in the Mnl20 clone. Furthermore, 5 out of 8 
Ncol clones analyzed lacked the Ugtl band (Figure 5B, Nco8, 
Ncol5, Ncol9, Nco20 and Nco27). These results confirmed the 
presence of a bi-allelic mutation in the clones Nco8, Ncol 5, 
Ncol 9, Nco20 and Nco27. The other clones analyzed showed 
reduced levels of Ugtl protein compared to the wt protein extract. 

We examined Ugt glucuronidation activity using the UGT-Glo 
Assay, by determining the percent consumption of the UGT 



multienzyme substrate, as described in "Materials and Methods". 
As controls, we used liver microsomes prepared from WT and 
Ugtl mutant mice [13]. Mutant mice have no Ugtl glucuronida- 
tion activity due to a targeted mutation in the Ugtl gene [13]. As 
expected, substrate consumption was significantly reduced in the 
Ugtl KO control (Figure 5C). The residual glucuronidation 
activity was performed by Ugt2, as the substrate used does not 
discriminate between Ugtl and Ugt2 activities. Since N-Muli cells 
express much lower amounts of Ugtl than mouse liver (Figure 5B, 
compare lanes 1 and 3, which correspond to 30 (0.g and 100 |Xg of 
total protein from WT liver and N-Muli cells, respectively), we 
performed the glucuronidation reaction using a higher amount of 
microsomes (10 u.g from cells and 1 |ig from liver). The amount of 
consumed substrate by N-Muli cells was much lower than in liver 
microsomes, confirming the lower expression levels of Ugtl in 
those cells. We observed that the biallelic mutant clones analyzed 
(Mnl20 and Nco20) had lower activity than N-Muli cells 
(Figure 5G). The activity of the monoallelic mutant clone analyzed 
(Nco 25) was intermediate between WT N-Muli cells and the 
biallelic clones. 
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Enzyme: 

TALEN: 



Uncut Ncol 



0 «V o<V 




B 



Ncol 



Ncol TALEN 




GFP 



- 



% of Uncut 100 so 



100 



<1 



% of Uncut 100 100 <1 8 



GGCATTCATCACACACTCTGGCTCCCATGGTATTTATGAAGGAATATGCAATGGAGTTC C57B1/6 

GGCATTCATCACACACTCTGGCTCCCATGGTATTTATGAAGGAATATGCAATGGAGTTC N-Muli 

GGCATTCATCACACACTCTGACTCCC GAAGGAATATGCAATGG AGTTC ( All ) 

GGCATTCATCACTCACTCTGGCTCCCATG TTATGAAGGAATATGCAATGG AGTTC ( A4 ) 

GGCATTCATCACACACTCTGGCTCCC TTATGAAGGAATATGCAATGG AGTTC ( A7 ) 

GGCATTCATCACACACTCTGGCTCCCAcGGTATTTATGAAGGAATATGCAATGGAGTTC ( +1 ) 

GGCATTCATCACACACTCTGGCTCCC TTATGAAGGAATATGCAATGGAGTTC ( A7 ) 

GGCATTCATCACACACTCTGGCTCCCAcGGTATTTATGAAGGAATATGCAATGGAGTTC ( +1 ) 

GGCATTCATCACACACTCTGGCTCCCAT TTATGAAGGAATATGCAATGGAGTTC (A5) 

GGCATTCATCACACACTCTGGCTCCC TGAAGGAATATGCAATGGAGTTC (A10) 

GGCATTCATCACACACTCTGGCTCCCAT TATGAAGGAATATGCAATGG AGTTC ( A6 ) 

GGCATTCATCACACACTCTGGCTC GTATTTATGAAGGAATATGCAATGGAGTTC ( A5 ) 

GGCATTCATCACACACTCTGGCTCCCATG TTTATGAAGGAGTATGCAATGGAGTTC (A3) 

GGCATTCATCACGCACTCTGGCTCCCATGtTATTTATGAAGGAATATGCAATGGAGTTC ( +1 ) 



Figure 3. Transient transfection of Ncol TALENs into N-Muli murine hepatoma cells. A) Determination of TALEN activity by RE-digestion of 
genomic PCR product. N-Muli cells were transfected with the Ncol TALEN pair together with a GFP plasmid. GFP positive cells were sorted, genomic 
DNA prepared and PCR amplified. PCR products were digested or not with Ncol, and the digestion product run in a 2% agarose gel. The arrow 
indicates the Ncol-resistant fraction of the PCR product to the Ncol-TALEN transfected cells (Lane 5). B) Confirmation of the presence of Ncol- 
resistant fragment. The agarose regions corresponding to the full-length undigested fragment of the gel shown in Panel A (GFP and Ncol-TALEN, 
lanes 4 and 5, respectively) were cut and DNA purified. A nested PCR reaction using internal primers was performed, and the PCR product was 
digested or not with Ncol (lanes 3, 6 and 2, 5, respectively). The percentage of uncut fragment is indicated bellow the gel. C) INDEL analysis of 
Ncol-TALEN N-Muli treated cells. The Ncol-resistant fragment of Panel A, Lane 5 was purified, cloned into a pUC19 vector and the obtained 
clones were sequenced. The left and right arms of the TALEN pair are indicated in green and red, respectively. The Ncol restriction enzyme site is 
indicated in blue. The deletions are indicated as black dashes. The type and length of the modification are indicated on the right of each sequence (A 
indicates deletion, and + base substitution, indicated by small orange caps). 
doi:1 0.1 371 /journal.pone.01 0481 6.g003 



We observed a good correlation between the results obtained 
with the reporter vector and the number of mutated clones 
obtained in vivo. In fact, the TALEN-pairs that showed to be 
highly efficient in the generation of INDELs in vivo, such as those 
targeting the Ncol and Mnll restriction sites, were also the most 
efficient ones in re-activating luciferase activity of the reporter 
vector (Table 1). Finally, the most efficient pair was the Ncol- 
TALEN, as determined by the Ncol-digestion assay of cell clones 
(Table 1), by the luciferase activity using the reporter vector 
(Figure 2D and Table 1), and by the number of clones lacking 
Ugtl protein TALENs treatment (Figure 5B and Table 1). 



Discussion 

The recent development of highly specific endonucleases 
opened the possibility to rapidly and efficiently modify the genome 
of cells, allowing the generation of cell lines having mutations in 
the gene of interest. These tools are of crucial importance to more 
deeply understand the mechanisms of disease and the metabolic 
pathways affected when specific genes are mutated. In the present 
work we generated specific TALEN pairs and tested their 
endonuclease efficiency in N-Muli murine liver-derived tissue 
culture cells, with the final outcome of creating mono- and bi- 
allelic mutations of the Ugtl gene. The rationale of using murine 
N-Muli cells instead of human cells to produce mutant cell lines 
resides in the future objective of performing in vivo editing of a 
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Table 1. Generation of Ugtl N-Muli mutant cells by TALENs. 





TALEN 
name 


Total N. of clones 
isolated 


Total N. of clones 
analyzed 


Total N. of RE- 
resistant clones 


TALEN efficiency (% of 
mutated cell clones) 


TALEN efficiency (obtained with the 
reporter CFP vector, relative to NlalV) 


NlalV 


10 


10 


0 


0 


1.0 


Ncol 


77 


30 


11 


36.6 


2.32 


BstXI 


12 


12 


2 


16.6 


1.03 


Mnll 


72 


30 


9 


30 


1.62 



The total number of clones isolated, analyzed and containing the RE site modified are indicated for all TALEN pairs used. The comparison between the TALEN pair 
efficiency based on the frequency of mutated cell clones obtained and that obtained with the reporter vector is shown. 
doi:1 0.1 371 /journal.pone.01 0481 6.t001 



Crigler-Najjar mouse model [13]. This mouse strain, which 
contains a one-base deletion in the exon 4 of the Ugtl gene 
resulting in a frameshift and premature stop codon, represents and 
excellent model to set up gene-editing using TALEN technology. 

We report a fast and efficient method to perform gene knockout 
in tissue culture cells, based on the use of TALEN pairs targeting 
restriction enzyme sites in the region of interest, as already 
performed in other systems [18-20]. The strategy of using 
restriction enzyme sites was necessary due to the presence of 
strain-specific allelic variants and pseudogenes [15,17], a condition 
that impeded the use of the Surveyor assay to detect INDELs. In 
addition, we presented a luciferase vector designed to rapidly and 
accurately determine the efficiency of the generated TALEN. 

The efficiency of generating double strand breaks and INDELs 
by ZNF, TALEN and CRISPR/CAS9 systems is usually 
determined by the mismatch assay (also called Surveyor, CEL-1, 
or Nuclease S assays), a method based in the detection and 
digestion by an endonuclease of mismatches in heterodimers 
formed between PCR-amplified DNA fragments from control and 
treated samples [21]. However, performing this technique in 
animal models faces the problem of the presence of strain-specific 



polymorphisms. Most targeted mutations are performed in 129/ 
Sv-derived ES cells, but mutated alleles are then transferred by 
successive backcrossing into a more widely used strain such as 
C57B1/6, which is one of the best characterized inbred strains and 
the reference strain for the mouse genome sequence database [22- 
24]. Consequendy, the wt gene of interest originates from a 
different genetic background compared to the mutant one. Thus, 
strain-specific allele polymorphisms may limit the use of the 
Surveyor assay to determine the efficiency of the gene editing/ 
gene KO, if strain-specific alleles are present. In addition, the 
presence of pseudogenes may complicate the analysis even more. 
We have shown here that these issues could be avoided by 
targeting restriction enzyme sites in the region of interest. To 
disrupt the Ugtl gene in N-Muli cells, we constructed TALEN 
pairs targeting four different restriction sites present in Exon 4. 
Despite the use of the same parameters to design the different 
TALEN pairs, our luciferase reporter assay showed that the 
cutting efficiency among enzymes varied substantially, suggesting 
the need of testing diverse TALEN-pair combinations in order to 
find the most efficient one. The causes for the variability are not 
clear, but may be related to the array length and composition [1 1]. 
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Figure 4. Analysis of N-Muli Ugtl -mutant cell clones. A) and B) Schemes of the PCR fragments showing the position of the PCR primers and 
restriction sites. The position of the Mnll and Ncol sites are indicated (Panels A and B, respectively), as well as the expected length of the digestion 
products. C) and D) Restriction enzyme assay to detect INDELs. DNA from antibiotic-resistant clones (Panels C and D, from Mnll and Ncol 
TALEN-transfected cells, respectively) was prepared, PCR amplified, digested or not with Mnll and Ncol and the products run in a 1 % agarose gel. The 
expected and obtained gel bands are indicated with arrows. N-Muli corresponds to control DNA, Mnl4, Mnl20, Nco25 and Nco20 to antibiotic- 
resistant clones that were transfected with the respective TALEN pairs. 
doi:1 0.1 371 /journal.pone.01 0481 6.g004 
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Mnll TALEN Left arm Mnll TALEN Right arm 

C57B1/6 WT TTTGGAAAATGCCCTTAAAACTGTCATCAACAACAAGAGgtaggtaatagaaaataaaaaacaaataaaaccaagaaccc 

N-Muli WT TTTGGAAAATGCCCTTAAAACTGTCATCAACAACAAGAGgtaggtaatagaaaataaaaaacaaataaaaccaagaaccc 

Mnll CI 4 TTTGGAAAATGCCCTTAAAACTGTCATCAACACCAAGG taatagaaaataaaaaacaaataaaaccaagaaccc 

Mnll CI 7 TTTGGAAAATGCCCTTAAAACTGTCATCAA 

Mnll CI 15 TTTGGAAAATGCCCTTAAAACTGTCATCAACA Ggtaggtaatagaaaataaaaaacaaataaaaccaagaaccc 

Mnll CI 20 TTTGGAAAATGCCCTTAAAACTGTCATCAACAACA---Ggtaggtaatagaaaataaaaaacaaataaaaccaagaaccc 

Mnll CI 20 

Mnll CI 24 TTTGGAAAATGCCCTTAAAACTGTCATCAACAACA-G-Ggtaggtaatagaaaataaaaaacaaataaaaccaagaaccc 

Mnll CI 29 TTTGGAAAATGCCCTTAAAACTGTCATCAA Aggtaggtaatagaaaataaaaaacaaataaaaccaagaaccc 



(A6) 

(A64) 

(A6) 

(A3) 

(A353) 

(A2) 

(A7) 



B 



Ncol 

C57B1/6 WT 
N-Muli WT 
Ncol CI 8 
Ncol CI 8 
Ncol CI 11 
Ncol CI IS 
Ncol CI 15 
Ncol CI 19 
Ncol CI 19 
Ncol CI 20 
Ncol CI 20 
Ncol CI 25 
Ncol CI 26 
Ncol CI 27 
Ncol CI 27 



Mouse 
Liver 



TALEN Left arm Nco I TALEN Right arm 

ACTCGGGCATTCATCACACACTCTGGCTCCCATGGTATTTATGAAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 
ACTCGGGCATTCATCACACACTCTGGCTCCCATGGTATTTATGAAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 

ACTCGGGCATTCATCACACACTCTGGCT TATGAAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 

ACTCGGGCATTCATCAC ACACTCTGGCT 

ACTCGGGCATTCATCACACACTCTGGCTCCCATG TGAAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 

ACTCGGGCATTCA - ATGGAGTTCCGATGGTGATGATGCC 



ACTCGGGCATTCATCAC ACACTCTGGCTCCCATG- 
ACTCGGGCATTCATCACACACTCTGGCTCCCA- - - 



AAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 

-ATGAAGGAATATGCAATGGAGTTCTGATGGTGATGATGCC 
CC 



ACTCGGGCATTCATCAC AC ACTCTGGCTCCC - - - 
ACTCGGGCATTCATCACACACTCTGGCTCCCAT- 
ACTCGGGCATTCATCAC ACACTCTGGCTCCC - - - 
ACTCGGGCATTCATCAC ACACTCTGGCTCCC - - - 
ACCCGGGCATTCATCACACACTCTGGCTCCCAT- 



- --TGAAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 
-TATGAAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 
-TATGAAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 
-TATGAAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 
-TATGAAGGAATATGCAATGGAGTTCCGATGGTGATGATGCC 



Mnl I clones 



Nco I clones 



4 7 15 20 24 29 8 11 15 19 20 25 26 27 



Ugtl 








(All) 

(A134) 

(A7) 

(A42) 
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Figure 5. Molecular analysis of cell clones. A) INDEL analysis of Mnll-TALEN and Ncol-TALEN N-Muli treated cells. The left and right arms of the 
TALEN pair are indicated in green and red, respectively. The restriction enzyme sites are indicated in blue. Deletions are indicated as dashes. The 
length of the deletions is indicated on the right of each sequence. B) Western blot analysis of the isolated cell clones. Proteins were prepared 
from the isolated cell clones and analyzed by Western blot using an anti-Ugtl antibody. Lanes 1 and 2 correspond to protein extracts prepared from 
WT and mutant livers [1 3], Lane 3 corresponds to untreated N-Muli cells. The arrow indicates the band corresponding to Ugtl . p-tubulin was used to 
normalize for protein loading. C) Ugtl glucuronidation activity of the isolated cell clones. Ugtl glucuronidation activity was determined by 
the UGT-Glo Assay as described in the Materials and Methods section. The activity determination was determined with microsomes prepared from wt 
and Ugt mutant mouse liver (1 ug), from wt N-Muli cells (10 |ig) and from the isolated cell clones (10 |ig, Mnl20, Nco20 and Nco25), as the percentage 
of luminescent substrate consumed. "+cntrl" indicates microsomes provided by the manufacturers (Promega) as positive control of the experiment. 
The experiment was repeated twice, and the mean of both experiments is shown. 
doi:1 0.1 371 /journal.pone.01 0481 6.g005 



Though, the two TALEN pairs that most efficiently rescued 
Luciferase activity of the reporter vector (Mnll and Ncol TALEN 
pairs) were also the most efficient ones in generating mutations in 
the Ugtl exon 4 in the genome context of N-Muli cells, obtaining 
a success rate of about 30% and 37% of mutated clones with those 
enzymes, respectively. Therefore, we reported a very good 
correlation between luciferase activity of the reporter vector and 
the number of clones having mutations in the exon 4 of the Ugt 1 
gene for the four TALEN pairs tested. Several cell clones with bi- 
allelic modification of the Ugtl gene were devoid of Ugt 1 -specific 
band in the Western blot analysis, confirming the drastic nature of 
the generated mutations. In the case of the Mnll site, it coincides 
with the 5' splicing site of Ugtl exon 4. Consequendy, it is 
expected that mutations in this site may affect pre-mRNA 
processing [25]. However, semi-quantitative RT-PCR failed to 
detect important differences in Ugtl mRNA levels (data not 
shown). More experiments should be performed to determine 
whether cryptic sites nearby are activated. In the case of the Ncol 
site, most of the analyzed mutations disrupted the natural reading 



frame, confirming the WB results where most clones contained bi- 
allelic mutations. Truncated proteins lacking the C-terminal 
portion lose the transmembrane-docking domain, resulting in an 
inactive protein [4,13,26-28]. In addition, despite of the fact that 
the TALEN-encoding plasmids did not confer any selective 
advantage to the cells (the selection markers were encoded by 
different plasmids), we cannot not exclude the eventual integration 
of the TALEN-encoding plasmids in the genome and further 
experiments should be performed to rule out this possibility. 

The study of glucuronosyl transferase enzyme activity was 
originally performed in microsomes from patients' liver [4,28,29] 
and is now normally carried out in tissue culture cell lines, such as 
COS-1, COS-7, fibroblasts or HEK293 ceU lines, stably or 
transiently transfected with plasmids expressing the specific 
cDNAs [28,30-33]. These cell lines are used to study substrate 
specificity [32], the requirement of specific co-factors and the 
potency of inhibitors [34] , and to correlate enzyme activity to the 
presence of specific mutations in the gene [33]. However, one 
limitation of this approach is that they do not generally derive 
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from the original tissue expressing the enzymatic activity. In that 
respect, the cell lines generated in the present work lack Ugtl 
glucuronidation activity and could be a very useful tool to perform 
the studies mentioned above and other pharmacological aspects in 
a more natural context, by adding back to the mutant cell clones 
specific variants of the Ugtl gene. 

To summarize, the mutant liver-derived cell lines generated in 
the present work represent a step forward to study biochemical 
aspects of Ugtl enzyme activity in a more natural context. We 
were able to generate efficient TALEN pairs and to determine the 
most efficient TALEN pair despite of the presence of strain- 
specific polymorphisms and pseudogenes. This step, needed to 
target the Ugtl gene, is a fundamental one towards setting up the 
gene editing technology, in order to find alternative therapies to 
treat the Crigler-Najjar syndrome type I. 

Supporting Information 

Figure SI Surveyor analysis of plasmid and genomic DNA 
sequences. PCR products of WT and MUT plasmids, which 
contain the murine Ugtl Exon 4 and flanking sequences differing 
by one nucleotide [13], were treated or not with Surveyor (lanes 2- 
3 and 12-13, respectively). Cleavage of the heterodimers (Lane 4) 
resulted in the expected DNA fragments of 192 and 324 bp. PCR 
products from WT and Ugtl -MUT mice were mixed and treated 
or not with Surveyor (Lanes 7 and 11, respectively). 
(TIF) 

Figure S2 Sequence of the WT exon 4 and mutation analysis by 
Surveyor. Alignment of reference GenBank sequence (C57B1/6 
strain) of WT exon 4 (capital letters) and flanking intron sequences 
(lower case) with that of 129/SvJ clones. DNA from 129/SvJ WT 
ES cells was PCR amplified, the PCR product cloned into pUC19 
and sequenced. The differences between the GeneBank and 129/ 
SyJ sequence are indicated with the red boxes. The SNPs present 
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